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Abstract 

To explain the effect of electron density saturation at high phos- 
phorus concentrations the model of negatively charged phosphorus 
clusters was compared with experimental data. A number of nega- 
tively charged clusters incorporating a point defect and phosphorus 
atoms (DPi)-, (DPs)", (DPa)^-, (DP3)-, (DPg)^-, (DP4)-, and 
(DP4)^~ were investigated at temperatures of 920 and 1000 °C. It 
was established that for clusters incorporating more than one phos- 
phorus atom the calculated electron density reached a maximum value 
and then decreased monotonically and slowly with an increase in the 
total dopant concentration. If a cluster incorporates 3-4 phosphorus 
atoms, the calculated dependencies of carrier concentration agree well 
with the experimental data for both temperatures regardless of the 
negative charge state of the cluster (— or 2—). Moreover, for all dou- 
bly negatively charged clusters investigated the fitting parameter is 
independent of the temperature. It means that formation of doubly 
negatively charged clusters is most likely but it is difficult to choose 
between clusters incorporating two, three or four phosphorus atoms. 
From the good agreement of the calculated curves with the experi- 
ments it follows that the model based on the formation of negatively 
charged clusters can be used in simulation of high concentration phos- 
phorus diffusion. 

Keywords: clusters; diffusion; doping effects; phosphorus; silicon 

1 Introduction 

At the present time, low energy high fluence ion implantation is widely used 
for producing the active regions of modern integrated microcircuits. A com- 
bination of ion implantation having energies of about 10 keV with rapid ther- 
mal annealing allows one to create doped layers with a depth of about 100 
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nm m 12] • To calculate the distributions of electrically active dopant atoms 
in such small regions, models of clustering having a high degree of accuracy 
should be used. The implementation of the adequate models of clustering 
is very important because saturation of electron density is observed experi- 
mentally for high dopant concentrations (see, for example [3]). Thus, only a 
part of dopant atoms is electrically active, and we have difficulties in achiev- 
ing high carrier concentrations. Unfortunately, all up-to-date models cannot 
adequately explain this phenomenon. Really, it is commonly accepted that 
the model of Tsai et al. was the first to account for saturation of charge 
carriers at an increasing doping level 0. In the model of ^4J it is supposed 
that arsenic clustering occurs as a result of the reaction 

3As- + e- ^"""^ As^ As3, (1) 

where As"*" is the substitutionally dissolved arsenic atom participating in 
cluster formation; e~ is the electron. 

The main feature of the model of is the assumption that clustered As 
atoms are electrically active (positively charged) at annealing temperatures, 
while at room temperature they are neutral. Therefore, the electron concen- 
tration at room temperature ur is approximately equal to the concentration 
of substitutionally dissolved arsenic atoms C. The assumption that As clus- 
ters are capable of trapping the conduction electrons and becoming neutral 
during cooling allows one to explain the saturation of electron density when 
— >^ cx) in in] . Here C'^ = C + C^^ is the total concentration of dopant 
atoms; C^'~^ is the concentration of dopant atoms incorporated into clusters. 
In Fig. 1, the electron density at room temperature nn calculated according 
to the model of Tsai et al. |3] is compared with the experimental data ob- 
tained by Nobili et al. The value of the fitting parameter in the mass 
action law for reaction (P) was chosen to satisfy the experimentally observed 
value of electron density saturation Umax which was equal to 2.11x10^ ^m~^ 

m 

As can be seen from Fig. 1, the curve calculated disagrees with the exper- 
imental dependence of electron density on the total dopant concentration. In 
Fig. 1, the dependences n ^ riR calculated for the neutral clusters VAs2 and 
VAs4 are also presented. Here n = n{C^) is the electron concentration at the 
temperature of cluster formation. From the calculations presented in Fig. 1 
it follows that when neutral clusters are formed, the electron density increases 
monotonically and saturation is not observed. Thus, dependences of electron 
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Figure 1: Calculated electron concentrations vs. the total arsenic concentra- 
tion for different models of clustering: solid line - model of Tsai et al. (Ref. 
jl]), dash-dotted line - neutral VAs2 clusters , and dashed line - neutral VAs4 
clusters. The experimental data (circles) are taken from Nobili et al. jS] 
for diffusion at a temperature of 900 °C. Dotted line represents the value of 
saturation. 

density calculated for neutral clusters also disagree with experimental data. 

In the clustering phenomenon was investigated with the example of 
Ge doped by As. It was proposed that cluster formation occurred due to the 
reaction 

mA+ + (m-ZCT)e" + D^< — >(AD)^^^^ (2) 

where A+ is the substitutionally dissolved donor atom; is the intrinsic 
point defect in the neutral charge state; m is the number of substitutionally 
dissolved atoms participating in the reaction; (AD)^c'' is the cluster formed; 
zci is the charge of this cluster. 

It is supposed that a cluster has the same charge during annealing and 
at room temperature. Based on this assumption, the dependence between 
the total dopant concentration and ur was obtained in [7j from the mass 
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action law for reaction Q 



im — 2;ctJ/ ^T\m 

C- - = K(T) "« , ' ' r . (3) 

{m-zci)'"' ^ 

where K{T) is the constant of local equilibrium for reaction (j21). 

Two different cases of cluster formation were investigated analytically 
based on dependence l^: (i) formation of neutral clusters {zci=0) and (ii) 
formation of negatively charged clusters including one dopant atom and a 
lattice defect. It was found that for neutral clusters the dependence = 
ur^C'^) has the following form: 

lgn. = -^lg(mA') + i^lgC-. (4) 

It follows from expression Q that at high impurity concentrations the 

line with the angle — in the logarithmical coordinates is asymptotic for 
2m 

the dependence hr = nji{C^) [T. This analytical result agrees with the 
numerical solutions presented in Fig. 1. 

In the case of a negatively charged cluster which incorporates one impurity 
atom and a lattice defect Eq.Q yields 

1 I I T^/-n^\ U + \^Cl\) 

It follows from Eq.(jSI) that saturation of electron density occurs if the 
dopant concentration approaches infinity. The maximum electron concentra- 
tion URmax has the following value (Tj : 



max 



\zci\K{T)] 1 + . (6) 

Thus, the model of 7j was the first to account for saturation of charge 
carriers at an increasing doping level. 

In the paper by Velichko et al. [H], the clustering of phosphorus atoms 
was investigated numerically. Due to numerical solution, the limitations of 
[7j were removed. It was proposed that the formation of negatively charged 
phosphorus clusters occurred due to the reaction 

mP+ + mz5D'^^ + A;e~^ — ^(PD)''^', (7) 



4 



and that the charge of a phosphorus cluster does not change during coohng. 
Here P"*" is the substitutionally dissolved phosphorus atom; D is the point 
defect participating in clustering; (PD)^'-^' is the cluster formed; m, rrir), and 
k are the numbers of impurity atoms, defects, and electrons participating in 
the cluster formation, respectively; ro and rci are the charge states of the 
defect and cluster, respectively. 

The possibility of formation of negatively charged phosphorus clusters is 
indirectly confirmed in ^ jTUl • It follows from calculations US] that singly 
(VAss)" and doubly (VAs2)^~ negatively charged clusters can be formed at 
high arsenic concentrations. 

The reaction used in fSj can be generalized to take into account the pos- 
sibility of defect generation during clustering: 

mP+ + miD[i + ke-< — ^(PD)'^^' + msDa^, (8) 

where D^^ and are the defects providing cluster formation and being 
generated during clustering, respectively; nii and m2 are the numbers of the 
defects participating in reaction ((Tj); ri and r2 are the charge states of these 
defects, respectively. For example, can be the self-interstitial generated 
during clustering [TTj . 

The charge conservation law is valid for chemical reaction (jSJ. This con- 
servation law is of the form 

m + rriiZi - k = zci + m2Z2, (9) 

where zi and Z2 are the charges of defects D^"^ and Dg^, respectively. 

For reaction (jH)) and conservation law ©, the equations obtained in jS] 
to describe clustering are modified to the following form: 

C^^ = KCn X^~^' rm + k + ^2^2) (j^n^ (iq) 



X 



1 

2ni 



m 
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where 



^ mA-.(/,g.)""n<-^-'"' + '= + ^''"''(Cg-xr. 

^"(S)^' ^'^^ 

Here x is the electron density normahzed to the concentration of intrinsic 
charge carriers in a semiconductor during diffusion ni\ C^'" and are 
the concentration of clustered phosphorus atoms and total concentration of 
acceptors, respectively. The parameter K has a constant value depending on 
the temperature of diffusion. This value can be extracted from the best fit 
to the experimental data. The parameters and li^^ ^^^^ the constants for 
the local thermodynamic equilibrium in the reactions when neutral defects 

and D2 are converted into charge states ri and r2, respectively; C^^^ 
and C^^^ are the concentrations of defects and D2 , respectively. The 
quantities C^^^ and C^^^ represent the equihbrium concentrations of these 
defects in the neutral charge state. 

In |H|, the formation of singly negatively charged clusters incorporating 
one phosphorus atom or two phosphorus atoms was investigated. It was 
shown that in the first case the concentration of charge carriers reached sat- 
uration if the concentration of substitutionally dissolved dopant atoms and 
correspondingly, the total concentration of phosphorus atoms had increased. 
Thus, these calculations agree with predictions of 7J. In the second case, the 
electron density reached a maximal value and then monotonically decreased. 
Unfortunately, in this paper a detailed comparison of calculated distribu- 
tions of electron density with experimental data was not carried out. Nor 
other cluster species were considered. Therefore, the purpose of this study is 
to investigate systematically the formation of different types of phosphorus 
clusters and compare the calculated dependences of electron density with 
experimental data. 



2 Calculations 

The system of equations (UHl) , (HH) obtained in jH] and modified in this paper 
allows one to consider different types of negatively charged clusters. In the 
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present investigation the electron density ur ^ n = n{C'^), concentration 
of substitutionally dissolved phosphorus atoms C = C{C^), and the concen- 
tration of clustered phosphorus atoms C^'" = C^'^{C'^) depending on the 
total phosphorus concentration were calculated for the clusters (VPi)~, 
(VP2)-, (VPa)'-, (VP3)-, (VPs)'-, (VP4)-, and (VP4)'-. To calculate 
these functions, a numerical solution of the system of nonlinear equations 
firU]) . (fTT|) was obtained by Newton's method. 

For comparison with the experimental data of and ^31; the process- 
ing temperature was chosen to be 920 °C {rii = 5.684x10^ /im and 1000 
°C {rii = 8.846x10^ /im ~^), respectively. In the experiments carried out in 
[12], the gas composition containing 0.27 % POCI3 was used as a source. The 
diffusion was carried out in a low dislocation density, (111) oriented, boron 
doped, Czochralsky pulled silicon with a resistivity of 1 cm. Duration 
of diffusion was 7 minutes at a temperature of 920 °C. The distribution of 
carriers concentration was extracted using the incremental sheet resistance 
and Hall measurements, performed after thinning of the specimens by anodic 
oxidation and oxide stripping. Neutron activation analysis was used for mea- 
suring the total concentration of impurity atoms. In silicon wafers, (100) 
oriented, p-type of 10 Q cm resistivity, were heavily implanted with phospho- 
rus at a fluence of 1.0x10^'' cm~^ and an energy of 100 keV and then with a 
fluence of 5.0x10^^ cm^^ at 50 keV. The samples were furnace annealed at a 
constant temperature of 1000 °C for 15 min. Secondary neutral mass spec- 
troscopy (SNMS) was employed to measure the dopant concentration profiles 
of the samples after annealing. Carrier concentration and mobility profiles 
were determined by the accurate incremental sheet resistance and Hall mea- 
surements, performed after thinning of the specimens by anodic oxidation 
and oxide stripping. 

In Fig. 2, the electron density calculated for the formation of (DPi)~ 
clusters is presented. As can be seen from Fig. 2, the saturation of electron 
density occurs at extremely high values of the total impurity concentration 
and calculated dependence ur ^ n = n{C^) disagrees with the exper- 
imental data. The disagreement is much less but also takes place for the 
formation of the (DP2)~ clusters. 

The electron densities calculated for the clusters (DP2)^^, (DPs)^^, 
(DPs)^, and (DP4)~ are presented in Figs. 3a, 3b, 4a, and 4b, respectively. 
As can be seen from these figures, a good agreement is observed between the 
calculated curves and experimental data especially for the case of (DP4)~ 
cluster formation. As follows from the calculated curves, for the clusters in- 
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Figure 2: Calculated concentrations of substitutionally dissolved phospho- 
rus atoms (solid line), clustered phosphorus atoms (dash-dotted line), and 
electron density (dotted line) vs. the total dopant concentration when singly 
negatively charged clusters incorporating a point defect and one phosphorus 
atom are formed. The experimental data (circles) are taken from Masetti et 
al. [12] for diffusion at a temperature of 920 °C. 

corporating more than one phosphorus atom the electron density reaches a 
maximum value and then decreases monotonically and slowly with increasing 

To calculate the electron densities presented, the generalized system of 
equations (fTUjl . (|TT|l was transformed due to assignment of the cluster charge 
state, a number of defects participating in cluster formation (one in the cases 
under consideration), and a number of phosphorus atoms incorporated into 
the cluster. For example, the following equations describing the concentra- 
tion of phosphorus atoms incorporated into the (DP4)~ clusters and electron 
density corresponding to the formation of these clusters can be obtained from 
the system dTTHl . (ITTll : 

C^"^ = KCdX' C\ (14) 
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(b) 

Figure 3: Calculated concentrations of substitutionally dissolved phospho- 
rus atoms (solid line), clustered phosphorus atoms (dash-dotted line), and 
electron density (dotted line) vs. the total dopant concentration for the for- 
mation of doubly negatively charged clusters incorporating two (a) and three 
(b) phosphorus atoms. The experimental data are taken from Masetti et al. 
Hg for diffusion at 920 °C. 
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Figure 4: Calculated concentrations of substitutionally dissolved phospho- 
rus atoms (solid line), clustered phosphorus atoms (dash-dotted line), and 
electron density (dotted line) vs. the total dopant concentration for the for- 
mation of singly negatively charged clusters incorporating three (a) and four 
(b) phosphorus atoms. The experimental data are taken from Masetti et al. 
for diffusion at 920 °C. 



It follows from the calculated dependence n = n{C'^) that with an in- 
crease in the total phosphorus concentration the concentration of charge 
carriers n reaches a maximum value Umax =3.25x10^ fim~^ at Cj^ax ~ 
6.90x10^ fim~^ and then monotonically decreases. At maximum the con- 
centration of substitutionally dissolved phosphorus atoms Cmax is 3.98x10^ 
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/iin ^. The fitting parameter KCd used in this calculation was equal to 
1.9x10-3^ yum^ 

For the case of doubly negatively charged clusters (DP4)^^ a disagreement 
is observed between the calculated dependence n = n{C^) and the experi- 
mental data. It follows from the comparison of the calculated dependences 
n = n{C^) with the experimental data that the suggestion on the forma- 
tion of negatively charged clusters (DP2)^", (DPa)^", (DP3)", and (DP4)~ 
provide a good fit of calculated functions to experimentally measured distri- 
bution of carrier concentration (see Figs. 3, 4). Therefore, it is difficult to 
extract the exact charge state and the number of the atoms incorporated into 
the cluster from the experimental data presented but the best fit is observed 
for the (DP4)~ clusters. 

In Figs. 5a, 5b, 6a, and 6b the dependences n = n{C'^) calculated for the 
clusters (DP3)", (DP4)~, (DPs)^", and (DP4)^~, respectively, are compared 
with the experimental data of obtained for a temperature of 1000 °C. 
As can be see from the figures, all these clusters give a good fit to the exper- 
imental data, especially the (DPa)^" and (DP4)^" clusters. For the case of 
formation of the (DPa)^" clusters, the system (jlO|) . (jlip has the form 

C"^^ = K CdX' C^, (16) 
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(17) 



and the value of the fitting parameter KCd was chosen to be equal to 
3.7x10-27 nm^. 

For the case of (DP4)2~ formation the system (fTn|l . (fTT|l has the form 
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Figure 5: Calculated concentrations of substitutionally dissolved phospho- 
rus atoms (solid line), clustered phosphorus atoms (dash-dotted line), and 
electron density (dotted line) vs. the total dopant concentration for the for- 
mation of singly negatively charged clusters incorporating three (a) and four 
(b) phosphorus atoms. The experimental data are taken from Solmi and 
Nobili for diffusion at a temperature of 1000 °C. 
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Figure 6: Calculated concentrations of substitutionally dissolved phospho- 
rus atoms (solid line), clustered phosphorus atoms (dash-dotted line), and 
electron density (dotted line) vs. the total dopant concentration for the 
formation of doubly negatively charged clusters incorporating three (a) and 
four (b) phosphorus atoms. The experimental data are taken from Solmi and 
Nobili ^31 for diffusion at a temperature of 1000 °C. 

and the fitting parameter KCn was chosen to be equal to 1.65x10"^^ fim^. 

As for (DP2)^~ cluster formation, there is a little disagreement between 
the calculated electron density and experimental data. It is interesting to 
note that for the case of formation of doubly negatively charged clusters 
incorporating more than one phosphorus atom the values of parameters KCd 



13 



providing a best fit to experimental data are equal for both temperatures 
920 and 1000 °C. At the same time, the values of KCu calculated for singly 
negatively charged clusters are not equal for different temperatures. For 
example, if a formation of (DP4)~ is assumed, then the values of KCd are 
equal to l-QxlQ-^s and 1.25x10-^^ urcv^ for temperatures 920 and 1000 °C, 
respectively. It follows from the extracted values of KCr, that the formation 
of doubly negatively charged clusters is most likely but it is difficult to choose 
between the clusters incorporating two, three or four phosphorus atoms. 
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3 Conclusions 



Analysis of the models dealing with the dopant atom clustering in silicon 
was carried out. It is shown that the models of clustering developed do not 
explain the available experimental data, especially the effect of electron den- 
sity satm^ation at high dopant concentrations. However, the phenomenon 
of saturation can be explained within the assumption of negatively charged 
phosphorus complex formation. To investigate this model, the electron den- 
sity and concentration of dopant atoms incorporated into clusters are cal- 
culated for temperatures of 920 and 1000 °C as the functions of the total 
phosphorus concentration. A number of negatively charged clusters incor- 
porating a point defect and phosphorus atoms, (DPi)^, (DP2)^, (DP2)^^, 
(DPa)-, (DPg)^-, (DP4)-, and {BF^f-, are investigated. It is shown that 
in the case of (DPi)" formation the concentration of charge carriers reaches 
saturation if the total phosphorus concentration increases. However, the sat- 
uration occurs only for very high values of the total dopant concentration 
that disagrees with experimental data. In the cases of clusters incorporating 
more than one phosphorus atom, the electron density reaches a maximum 
value and then decreases monotonically and slowly with an increase in the 
total dopant concentration. If a cluster incorporates 3-4 phosphorus atoms, 
the calculated dependences of carrier concentration agree well with the ex- 
perimental data for both temperatures regardless of the value of the negative 
charge state of the cluster (— or 2—). On the other hand, for all doubly 
negatively charged clusters investigated the values of the fitting parameter 
are independent of the temperature. It means that formation of doubly neg- 
atively charged chistcrs is most likely but it is difficult to choose between 
clusters incorporating two, three or four phosphorus atoms. From a good 
agreement of the calculated curves with experimental data it follows that 
the model based on the formation of negatively charged clusters can be used 
in simulation of high concentration phosphorus diffusion. 
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